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o ; Abstract 

The time-dependent barrier passage of an anomalous system-reservoir coupling non-equilibrium 
open environment is studied where the heat bath is modulated by an external noise. The time- 
' dependent barrier passing probability is obtained analytically by solving the generalized Langevin 

equation. The escaping rate and effective transmission coefficient are calculated by using of the re- 
active flux method in the particular case of external 5-correlated noise modulated internal Ornstein- 
Uhlenbeck process. It is found that not all the cases of external noise modulating is harmful to 
the rate process. Sometimes it is even beneficial to the diffusion of the particle. There exists an 
optimal strength of external noise modulation for the particle to obtain a biggest probability to 
escape from the potential well to form a notable rate of effective net flux. 
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I. INTRODUCTION 



The diffusion model for chemical reactions has become ubiquitous in many areas of 
physics, chemistry and biology since its advancing by H. A. Kramers in terms of the theory 
of Brownian motion in phase space |lj . Several of its variants have been proposed for under- 
standing the nature of activated processes in classical (2 - 5 ] , semiclassical jfj], and quantum 
systems [7 10 ] . The majority of these treatments concern essentially with an equilibrium 
thermal bath at a finite temperature that simulates the reaction coordinates to cross the ac- 
tivation energy barrier and the inherent noise of the medium originates internally. Therefore 
these systems are actually classified as a thermodynamically closed system-plus-reservoir en- 
vironment in contrast to the systems directly driven by external noise (s) in non-equilibrium 



statistical mechanics 



111 ] . However, what will happen if the total system-reservoir environ- 



ment is really open or modulated by some external forces, is the barrier crossing dynamics 
still amenable to the present theoretical analysis ? These problems are of valuable consid- 
eration from a microscopic point of view. 

Recently, A type of system-reservoir coupling model simulating an activated rate process 



is studied where the heat bath is modulated by an external noise [12|, [3]. A non-equilibrium 
stationary state distribution which is reminiscent of the equilibrium Boltzmann distribution 
in calculation of rate is proposed for this open environment by analytically solving the gen- 
eralized Fokker-Planck equation and a generalized Kramers' rate is derived via the method 
of flux-over-population |14j. This brings a new broad perspective for the study of anomalous 
activated processes in non-equilibrium open systems. However, some crucial problems such 
as the barrier recrossing phenomenon which is always met in an escaping process have not 
been concerned in these works. The time-dependent dynamical barrier escaping process in 
non-equilibrium open systems is far from well studied. We noticed that these can be easily 



achieve d by the method of reactive flux 



15Ml71] which is developed from the transition state 



theory 18M20I]. Therefore, following this powerful and convenient way we present in this 
paper a careful study on the non-equilibrium activated barrier crossing process concerning 
intensively on its dynamical details. 

The paper is organized as follows: In Sec. II, the expression of the barrier passing 
probability is obtained by analytically solving the relevant generalized Langevin equation. 
In Sec. Ill, we give the time-dependent barrier escaping rate and transmission coefficient 
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derived by using of the reaction flux method. Sec. IV serves as a summary of our conclusion 
where some implicate applications of this study are also discussed. 



II. MODULATED SYSTEM-RESERVOIR ENVIRONMENT 



The physical scenario depicting the modulation of the bath by an external noise lives in 
many different kinds of situations in forming a non-equilibrium (or non-thermal) system- 



reservoir coupling environment 2jJ, |22j. For example, the simple unimolecular conversion 



process from A to B in an isomerization reaction is generally considered to be carried out in 
a photochemically active solvent under the influence of external fluctuating light intensity. 
In the reaction, fluctuations in the light intensity result in fluctuations in the polarization 
of the solvent molecules. Thus the effective reaction field around the reactant system gets 
modified. Given the required stationarity of this non-equilibrium open system is maintained, 
the dynamics of barrier crossing evolves amenable to the present theoretical analysis that 
follows. 

Let us begin our study from the Hamiltonian describing a system of particles with unit 
mass bilinearly coupled to a harmonic reservoir that is modulated by an external noise. 
Mathematically it reads 23] 

N i i i 

1 +H int , (1) 
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H=-p 2 



U(x) + 



2^i + -^Mi- C i x f 



where {x,p} and {qj,Pj} are the sets of coordinate and momentum variables of system and 
reservoir oscillators, respectively. CjX measures the interaction between the particle and 
reservoir, U(x) is the potential. H int = J2j=i ^jQj 6 ^) represents the modulating interaction 
on the reservoir results from the external noise e(t) which is assumed to be stationary 
and Gaussian with zero mean and decaying second order correlation function (e(t)e(t')) e = 
2Dtp(t — t'). Here (■ ■ • ) e implies the averaging over all the realizations of e(t) with strength 
D and ip(t) is a relevant memory kernel. 
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Eliminating the bath variables in the usual way 
equation (GLE) can be obtained as 



x + dt'j(t-t')x(t') + d x U(x) = ((t)+s(t), 
Jo 



25j, a generalized Langevin 



(2) 
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where ((t) = £\ c j{[<?j(0) — CjX }oj^cosUjt + pj(Q)ujSmojj t} is the internal fluctuation force 
generated from the system-reservoir coupling. While e(t) = - /„ (p(t-t')e(t')dt' with <^(t) = 

■ CjXjUjSinUjt is an additional fluctuation force resulted from the external noise e(t). The 
anomalous form of Eq.(j2J suggests that the system is under the combining government of 
two forcing functions. This will in no doubt lead to some novelty results. 

Before reaching the central point of this article, let us firstly digress a little bit about 
C(t) and e(t). Due to its particular origin, the statistical properties of are determined 
by the initial conditions of the system-reservoir coupling environment which is assumed to 
be equilibrium at t — when the external noise agency has not been switched on. That is 



(at) 
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) = and {((t)((t')) = k B T^{t — t') satisfying the fluctuation and dissipation theorem 



27J. However the statistical properties of e(t) may be enslaved to several aspects of 



factors such as the normal-mode density of the bath frequencies, the coupling of the system 
with the bath, the coupling of the bath with the external noise and the external noise itself. 
The very structure of e(t) suggests that this forcing function is different from a direct driving 
force acting on the system. Therefore after the external noise agency is switched on, the 
system can be regarded as being driven by an effective noise (= ((t) + e(t)) whose 
correlation is given by 

(0» = k BTl(t - t') + 2D [ dtMt -h) I 1 dt^ih - t 2 )<p(t - t 2 ), (3) 

Jo Jo 

along with ((£(£))) = 0, where ((■••)) means taking two averages independently. This 
relation is reminiscent of the familiar fluctuation-dissipation theorem. However, due to the 
appearance of the external noise intensity, it rather serves as a thermodynamic consistency 
condition instead. 

Due to the Gaussian property of the noises ((t) and e(i) and the linearity of the GLE, 
the joint probability density function W(x,v ,t; xo,vo) of the system oscillator must still be 
written in a Gaussian form 

W(x,v,t;x ,v ) = ^i^ e -MvnOA^WvW], (4 ) 
where y(t) is the vector [x — (x(t)),v — (v(t))} and A(t) is the matrix of second moments 
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with each component as 



A n {t) = (ix-(x(t))] 2 ), 
A 12 (t) = A 21 (t) = ([x-(x(t))][v-(v(t))}), 
A 22 (t) = (iv-(v(t)}f). 

The reduced distribution function can then be yielded by integrating Eq. (j4]) over v as 



(5a) 
(5b) 
(5c) 



W(x,t;x ,v ) 



1 



zexp 



2A n (t) 



(6) 



in which the average position {x(t)) can be obtained by Laplace solving the GLE. In the 
case of an inverse harmonic potential U(x) 



■\mJlx 1 , it reads 



(x(t)) 



H(t')dt' 



x + H(t)v 



(7) 



where H(t) namely the response function can be yieldec 



H(s) = [s + S7(s) — u%] with residue theorem 



3QJ. 



from inverse Laplace transforming 



III. BARRIER ESCAPING PROCESS 



For the activated barrier crossing process, one of the most crucial factors that we concern 
is the probability of passing over the saddle point (namely also the characteristic function) 
which can then be determined mathematically by integrating Eq. ([6]) over x from zero to 
infinity 



x(x ,v ;t) = J W(x,t;x ,v )dx, 
2 I 



(8) 



which will lead to a finite real number range from to 1 with 1 for reactive trajectories 
and for nonreactive ones. The escape rate of a particle, defined in the spirit of reactive 
flux method by assuming the initial conditions to be at the top of the barrier, can then be 
yielded from 

" oo 



kit) 



h 



dx Q / v W st (x , v )5(x - x b )x{x , v ; t)dv 



(9) 
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in the phase space, where W s t(xo,v ) generally should be an equilibrium Boltzmann distri- 
bution that depends on the initial position and velocity of the particle. However, in the 
external noise modulated system-reservoir coupling environment it should be replaced by a 
Boltzmann form stationary probability distribution W st (xo,v ) = iexp[— {^y- + w+sst-^, }] 



2D b 1 D b +V(oo) 

which has been proved in Ref. [13] as a steady-state solution of the relevant Fokker-Planck 
equation. In the formula aforesaid, Q is the partition function, U{x) = U^ — §fi&(<E — Xb) 2 is 
the renormalized linear potential near the barrier top with f2& an effective frequency and C4 
the barrier height. \P(oo) and Dp are two asymptotic constants to be calculated in the long 



time steady state (refer to Ref.[13| for detailed information). This stationary distribution 
for the non-equilibrium open system is not an equilibrium distribution but it plays the role 
of an equilibrium distribution of the closed system, which may, however, be recovered in the 
absence of the external noise. 

In general, the above total rate k(t) in Eq.(j9]) can be viewed as a generalized TST rate 
k TST = -lj-e~ Ub ^ Db+ " s ' ( - 00 ^ (where [Db + \I/(oo)]//cb in the exponential factor defines a new 
effective temperature characteristic of the steady state of the non-equilibrium open system) 
multiplied by a factor between and 1 which describes the possibility of a particle already 
escaped from the metastable well to recross the barrier. By substituting Eqs.(j5]), © and 
© into Eq.® we obtain k{t) = n(t)k TST with 



D b H\t) 

acting as an effective transmission factor which leads immediately to the Kramers formula 

nn 

for the rate constant [ll. Il4j in the absence of the external noise. As expected, both n(t) and 
k(t) are the functions of the external noise strength D and the coupling of noise to the bath 
modes. The varying of K(t) provides an isolated inspection on the dynamical corrections of 
k(t) to the TST rate. 

In the following calculations we rescale all the parameters so that dimensionless quantities 
such as ksT = 1.0 are used. We consider a particular case as an example where the external 
noise is 5 correlated and the internal is an Ornstein-Uhlenbeck (OU) process, i.e., 

(e(t)e(t')} e = 2D5(t-t') } (11a) 
(C(i)C(O) = k B T^e-^\ (lib) 

T 

both symmetric with respect to the time argument and assumed to be uncorrelated with 



each other. Where D and g are effective friction constants and r _1 is a cutoff frequency of 
the system oscillator. It should be noted that for r — > 0, the internal noise shown above 
becomes also 5-correlated. Asymptotic parameters contained in the stationary probability 
distribution W st (x ,v ) are determined to be fij = oj b , D b = k B T and ^(oo) = in the 
particular case that we considered here. 



1.0 - 



0.8 



> 0.6 

o 

0.2 
0.0 



"9=3 
9=5 
9=7 
9=9 



1.0 

0.8 

^-,0.6 
-t— > 

2o.4 
0.2 
0.0 



"9=3 
g=A 
9=5 
9=6 



50 



100 



150 



200 



50 



100 



150 



200 



FIG. 1: (Color online) Instantaneous values of x( x o,vo;t) and n(t) for various strengths of purely 
internal noise (D = 0). Parameters used are r = 3.5, uj^ = 1.2 and xq = —1.0, vq = 0.5 respectively. 
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FIG. 2: Stationary values of Xat and function of the strength of internal noise g where 

identical parameters are used as those in Fig{TJ 

Firstly, we present a detailed investigation on the rate process which is purely determined 
by the internal noise. FigJT] gives the instantaneous values of the barrier passing probability 
x{xo,v ;t) and transmission coefficient n(t) for various strengths of internal noise g where 
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all the system parameters are set dimensionless and particles are assumed to start from 
Xo = —1.0 position with initial velocity vq = 0.5. From which we can see that both x{ x o, v o] t) 
and n(t) evolve asymptotically to a stationary value in the long time limit. Meanwhile the 
stationary values of them decrease monotonously as the increase of the strength of internal 
noise. This can also be witnessed in Figf2]where the stationary values Xst and K st are plotted 
as a function of the strength of internal noise g. This is a trivial phenomenon which can 
be understood in ease for a rate process. Because a hard dissipative environment is in no 
doubt harmful to the diffusing process. Not only will it present a frictional resistance force 
but also can it increase the amplitude of fluctuation. Therefore a smaller and smaller net 
flux is expected as the increasing of the strength of the internal noise. Particles which have 
already escaped from the potential well will have a big probability to recross the barrier. 
Correspondingly related is a small n st as is shown in Fig{TJ 
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FIG. 3: (Color online) Instantaneous values of xi x 0i v o> t) an d n(t) for various strengths of external 
noise. Parameters used are r = 3.5, w& = 1.2, g = 3.5 and xq = —1.0, vq = 0.5 for each case. 

In further, let us turn to the most important point of our study, i.e., what will happen if 
the rate process is additionally modulated by an external noise? In order to obtain a explicit 
elucidation of this problem we present in Figj3] again the instantaneous values of x{ x o, v o'i t) 
and n(t) at various strengths of external noise where identical system parameters are used 
as those in FigJT] except for D = 0, 2, 3, 4, 8 and 18 respectively. The stationary values of 
them (xst and /t st ) are also exhibited as a function of the external noise strength D in FigJH 
in the meantime. From which we can see a non-trivial phenomenon that both Xst and K st 
varies non- monotonously as the increasing of the strength of external noise. That is to say 
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FIG. 4: Stationary values of Xst and K s £ clS cl function of the strength of external noise where 
identical parameters are used as those in Figj3j 

not all the cases of external noise modulating is harmful to the rate process. Sometimes 
it is even beneficial to the diffusing of the particle. For it can result in a bigger Xst and 
/Cat comparing to the pure internal case. This reveals, in the combining effect of internal 
and external noises, the particle will be able to get a biggest probability to escape from the 
potential well. Therefore a notable rate of net flux is expected. In further, we can infer 
from it that there lives an optimal strength of external noise modulation for the particle 
to obtain a biggest barrier escaping probability. In the particular case that is considered 
here this optimal value is about D = 3.4 (dimensionless) as can be calculated following the 
procedures in the context foreknown. 

IV. SUMMARY AND DISCUSSION 

In summary, we have studied in this paper the time- dependent barrier passage of an 
external noise modulated system-reservoir coupling environment. We proposed analytically 
the barrier passing probability and the rate of barrier escaping. The influence of the exter- 
nal noise on the dynamical barrier escaping process is investigated carefully by calculating 
the effective transmission coefficient of an external (5-correlated noise modulated internal 
Ornstein-Uhlenbeck process for an example. The main conclusions of our study is that not 
all the cases of external noise modulating is harmful to the rate process. Sometimes it is 
even beneficial to the diffusion of the particle. There exists an optimal strength of external 
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noise modulation for the particle to obtain a biggest probability to escape from the potential 
well to form a notable rate of effective net flux. 

In applications and industrial processing, the creation of a typical non-equilibrium open 
situation by modulating a bath with the help of an external noise is not an uncommon 
phenomenon. The external agency generating noise does work on the bath by stirring, 
pumping, agitating, etc., to which the system dissipates internally. However, provided the 
long-time limit of the moments for the stochastic processes pertaining to the external and 
internal noises characterized by arbitrary decaying correlation functions exist, the expression 
for the effective transmission coefficient of barrier crossing rate for the open system we 
derive here is fairly general. Therefore, we believe that these considerations are likely to 
be important in other related issues in non-equilibrium open systems and may serve as a 
basis for studying processes occurring within irreversibly driven environments 
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for thermal ratchet problems 33]. The externally generated non-equilibrium fluctuations 
can bias the Brownian motion of a particle in an anisotropic medium and may also be used 
for designing molecular motors and pumps. 
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